Thermoelectric efficiency in the space-charge-limited transport regime in 

semiconductors 
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The thermoelectric efficiency of semiconductors is usually considered in the ohmic electronic 
transport regime, which is achieved through high doping. Here we consider the opposite regime of 
low doping where the current- voltage characteristics are nonlinear and dominated by space-charge- 
limited transport. We show that in this regime, the thermoelectric efficiency can be described by a 
single figure of merit, in analogy with the ohmic case. Efficiencies for bulk, thin film, and nanowire 
materials are discussed, and it is proposed that nanowires are the most promising to take advantage 
of space-charge-limited transport for thermoelectrics. 



I. INTRODUCTION 



Most implementations of thermoelectric devices uti- 
lize highly doped semiconductors with linear (ohmic) 
current-voltage characteristics. In such materials, the 
thermoelectric efficiency is determined by the figure of 
merit zT — S 2 aT/ k where S is the Seebeck coefficient, a 
the electrical conductivity, T the temperature, and k the 
thermal conductivity. The zT factor has governed much 
of the research in developing approaches to improve ther- 
moelectric efficiency; however, the intimate relationship 
between transport quantities places constraints on the 
possible paths towards efficiency improvements. Thus, 
new approaches where the thermoelectric efficiency does 
not depend on zT would open new routes for scientific 
exploration. 

An example of such an approach is thermionic cool- 
ing and power generation, which was originally con- 
sidered for injection into vacuum]!], and then into 
semiconductors ;2j. These thermionic devices are often 
considered different from thermoelectric devices due to 
the presence of an injection barrier and because the 
channel length is shorter than the electronic mean-free 
path, giving ballistic electron transport [3|-|5[- However, 
as the channel length increases beyond the mean-free 
path, thermionic devices become equivalent to thermo- 
electric devices [1, 0], which has led to new proposals to 
engineer multilayer semiconductor materials y|. 

In this manuscript, we explore an alternative approach 
that does not rely on ohmic transport or thermionic emis- 
sion, but instead exploits space-charge-limited (SCL) 
electronic transport in a semiconductor. We derive an 
expression for thermoelectric efficiency in this regime, 
and show that it depends on a new dimensionless fig- 
ure of merit that replaces zT . Furthermore, we apply 
the theory to bulk, thin film, and nanowire materials, 
and conclude that nanowires are the most promising to 
harness SCL transport for thermoelectrics. 

The term " space-charge- limited" originated from the 
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FIG. 1: Panels (a) and (b) show band diagrams for a metal- 
semiconductor-metal system when the semiconductor has low 
doping. Figure (a) is at equilibrium while figure (b) is un- 
der an applied voltage. Solid lines are the valence and con- 
duction band edges and dashed lines are the Fermi levels, 
(c) Power loss due to Joule heating plotted for ohmic trans- 
port and for SCL transport in dimensionless units. The inset 
shows the current-voltage characteristics for these two trans- 
port regimes. 



consideration of charge injection into a region where 
electronic transport is ballistic. In that case, space- 
charge effects are detrimental because they create an ad- 
ditional barrier for injection compared with the reference 
space-charge-free system[9|]. This situation is relevant to 
thermionic devices where space-charge effects have been 
shown to reduce efficiency [l|. The issue of space-charge 
is different if one considers semiconductor-based thermo- 
electric devices. There, the reference system is a highly- 
doped semiconductor with two ohmic contacts. In that 



2 



case, space-charge effects are minimaf because the high 
doping immediately screens the injected charge. How- 
ever, if the doping in the semiconductor is lowered, space 
charge effects become dominant, and the current becomes 
larger than the ohmic current because electronic trans- 
port is no longer determined by the free-carrier relaxation 
time, but by the shorter carrier transit timeQ. This situ- 
ation, which has been observed experimentally in a broad 
range of materials, is the one that we consider in this 
manuscript. 

Figure 1 shows the system under consideration: a low- 
doped semiconductor between two ohmic contacts. Band 
diagrams are shown in Figs la,b for the case relevant 
to electron injection: the metal Fermi level contacts the 
semiconductor in the conduction band creating a large 
density of carriers in the near-interface region. Band- 
bending away from the metal/semiconductor interface 
creates a barrier for electron injection at zero bias. Un- 
der bias, a maximum in the potential is created near 
the injec ting electrode, defining the so-called "virtual 
cathode" . Modeling has demonstrated that in the 
presence of this virtual cathode, SCL transport domi- 
nates as soon as the applied bias V exceeds a few fcr[TTj]. 
This leads to a current-voltage relationship J oc V 2 as 
discussed further below. This relationship differs from 
the exponential J — V characteristics usually considered 
for thermionic systems [lj]; the transitions between the 
thermionic and SCL regimes and the role of Schottky 
barriers has been studied in 

detailtHEDl, 

demonstrating 

that SCL dominates even when small Schottky barriers 
are present. 



as 



v(x) = (x/L) 3/2 V. 



(2) 



This nonlinear spatial dependence of the potential leads 
to the current-voltage relationship 



j=ev 2 . 



(3) 



Thus, in the SCL regime, the current is no longer linear in 
V, but is quadratic instead. The proportionality constant 
£ depends on the dimensionality of the material, and to 
some extent on the shape of the electrodes, as will be 
discussed further below. 

The quadratic form of the J — V behavior implies that 
P e i = J 3 / 2 /£; the interest in this transport regime for 
thermoelectrics is illustrated in Fig. lc, where P e i for 
the ohmic and SCL regimes are compared. There, P e i 
in the SCL regime is below that of the ohmic regime 
at sufficiently high voltage, suggesting a potentially new 
regime of operation for thermoelectrics. (The cross-over 
between ohmic and SCL transport occurs at V c = a/£ 2 L, 
as shown in the Fig. lc inset. This sets a condition for 
when ohmic contributions should be negligible compared 
to SCL transport.) 

To derive the thermoelectric efficiency in this regime, 
we utilize the approach described in Ref. [14|. From the 
above equations we obtain for the local efficiency 



drj[x) 
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T S + kVT/TJ 
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— Vr (x,T), (4) 



II. THERMOELECTRIC EFFICIENCY 

To be specific, we consider the thermoelectric power 
generation efficiency of the system of Fig. 1, but the 
approach should apply to cooling as well. The efficiency 
of a thermoelectric material for power generation is given 
by0 



where we used dP e i = Jdv(x), and where the prefac- 
tor dT/T indicates that the efficiency is limited by the 
Carnot efficiency. The temperature T depends on posi- 
tion through the heat equation 
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dx 
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V = 



Pth ^ = J Jl h S(x)dT - J J L Vv(x)dx 
Q JT h Sh + K/iVTh 



(1) 



where P t h is the thermoelectric power generated, P e i is 
the electrical power dissipated, and Q is the thermal 
power supplied. Here T c and Th are the cold and hot 
side temperatures and v(x) is the local potential. 

Usually, one considers a highly doped semiconductor of 
length L between ohmic contacts. For a bias voltage V, 
there is a linear drop of the potential along the channel 
due to the free-carrier screening; this leads to J = aV/L, 
giving P e i = J 2 L/a and the dependence of the efficiency 
on zTjljj. For SCL transport, the transport equations 
need to be solved self-consistently with the Poisson equa- 
tion due to the unscreened space-charge [15j; in the sim- 
plest model with diffusive electronic transport [l5j. this 
gives an electrostatic potential that depends on position 



The first term on the right hand side is the usual Thom- 
son effect, while the second term is the local Joule heating 
in the material dP e i . This equation applies for SCL trans- 
port because both the electronic and phononic transport 
are diffusive, and thus both the electronic and phononic 
thermal conductivities are well defined (this has been 
demonstrated for SCL transport in nanowires [l7j ) . 
The total efficiency of a segment is obtained from[l8| 



rj = 1 — exp 



T 



dT 



(6) 



Noticing that d In (ST + kVT/ J) / dT = [S + 
d(KVT/J)/dT]/[ST + kVT/J] and assuming that 
the coefficients are independent of temperature and 
position we arrive at 
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S C T C + kVTJJ 
S h T h + KVT h /J' 



(7) 



Thus, the efficiency depends only on quantities evaluated 
at the cold and hot sides[14]. To obtain these quantities 
we integrate the heat equation to get 



kVT c /J 
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which leads to 
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where 



J 



1/2 



and 



5 2 r (AT)' 



(11) 



The variable y is a dimensionless quantity equal to the 
ratio between the generated voltage and the maximum 
thermoelectric voltage that can be developed across the 
material. gT is a dimensionless parameter that replaces 
the zT factor from conventional ohmic materials. (Note 
that the average temperature T does not directly appear 
in gT. However, we use the notation T as a reminder 
that, as a first approximation, the quantities that enter 
gT should be evaluated at the average temperature across 
the material.) 

Figure 2 shows the dependence of the efficiency on 
y from Eq. ([5]), indicating that positive efficiency is 
achieved in the range < y < 1, with a maximum at- 
tained at an intermediate value. This behavior is similar 
to that obtained for conventional ohmic losses but 
with a different functional dependence. 

The maximum efficiency is obtained by maximizing r\ 
with respect to y to obtain 
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(12) 



This equation provides the optimal y (gT, AT/Th) that 
maximizes the efficiency: 



V 



1 A 1 



(13) 



where 5 = + 2V2%/TT2i] 1 and e = 

^A^1-|A^. At this optimal y the efficiency is 




FIG. 2: Thermoelectric efficiency as a function of variable y 
for three values of the gT factor. 



(10) given by 
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(14) 



Thus, the efficiency only depends on gT and AT/Th, 
much like ohmic materials where the dependence is on 
zT and AT/T h . 

Figure 3a shows the calculated efficiency as a function 
of gT for AT/Th = 1/2, indicating low efficiency at small 
gT and efficiencies approaching the Carnot efficiency at 
large gT. The two limiting cases of small and large gT 
can be obtained as 



AT 



" 3 (l loTf) {n9 T 
(ff)%T small ^gf 
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_\ -1/3' 



large ^gT 



(15) 



thus providing simple analytical expressions for these two 
regimes. The small (AT/Th) gT expression provides a 
good estimate of the full curve up to about 1% efficiency, 
while the large (AT/Th) gT approximation is valid for 
V £ 20%- 

The inset in Fig. 3a shows the efficiency for a tradi- 
tional ohmic material as a function of zT for AT/Th = 
1/2, where zT = 1 gives about 10% efficiency. Since this 
is the current state-of-the-art in thermoelectric materi- 
als, we can use it as a comparison with the SCL regime; 
in that case one would need gT ss 5 to achieve the same 
efficiency at this given AT/Th- However, the efficiency 
depends significantly on AT/Th, as shown in Fig. 3b for 
values of gT as high as 10. The curve for gT = 1 is well 
approximated by the small (AT/T h ) gT limit of Eq. flT5J| , 
indicating a quadratic dependence on AT/Th- 

The full expression for gT depends on the parameter 
£, which in turn depends on the dimensionality of the 
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FIG. 3: (a) Thermoelectric efficiency as a function of the gT 
factor. Dashed lines are the low and high gT behaviors from 
Eq. (|15|) . Inset is the efficiency as a function of zT. (b) 
Thermoelectric efficiency as a function of AT/T^. 



material under consideration; for bulk[l5j], thin film [I? 
and nanowire fl9| materials, it is given by 



|p bulk 

£ j^ thin film, t < L 
£o WT nan °wire, R <C L 



(16) 



where e is the permittivity, /i is the mobility, and £ is 
a numerical constant that depends on the shape of the 
electrodes, equal to 1 for planar contacts t is the 
thickness of the thin film, while R is the nanowire radius 
(see Fig. 4 for illustrations of the different geometries.) 
The expressions for the thin film and nanowire cases are 
applicable when i/L<l and R/L -c 1; for larger values 
oft/L and R/L,£ crosses over to the bulk expression p^j. 

The gT factor also depends on the material parame- 
ters S, k, /i, and s. Like traditional ohmic materials, 
high gT requires large S, low k and large /i; but it also 



needs large e, a criterion not usually required for high zT. 
This arises because large e serves to screen the injected 
charge, and reduces the repulsive Coulomb interaction 
that opposes charge injection. This is beneficial because 
dielectric screening maintains the SCL transport regime, 
in contrast to free-carrier screening that would make the 
current ohmic. 

Two advantages of the SCL regime are that the maxi- 
mum values of S and n can be exploited. Indeed, ohmic 
materials are usually operated at high doping where 
both S and fj, are reduced from their maximum values. 
But at low doping S attains a maximum at a value (20j 
S ~ E g /2eT while \i saturates to its intrinsic value \i int - 
We can thus write the maximum value of gT as 



2eT 



(17) 



where a is a geometry-dependent factor that can be ob- 
tained from Eq. (fl"6|). We note that the gT factor de- 
pends linearly on /i, just like an ohmic material where 
zT a ~ fj,- however, the difference is that the mobility 
in zT is the high-doping mobility, which is smaller than 

Mint ■ 



III. PRACTICAL ESTIMATES 

To evaluate the practicality of the SCL approach for 
thermoelectric power generation, we apply the theory to 
the bulk, thin film, and nanowire geometries for T^ = 600 
K and AT = 300K. We consider GaAs at low doping as 
an example material since SCL has been observed in this 
material|21], and because material properties are readily 
available [22]: at T = 450K and low doping, E g =1.35 



eV, e 



/' 



3150 cm 2 /VS, and 



30 W/mK. 



Figure 4a shows the calculated efficiency for a bulk device 
as a function of the channel length L. The efficiency is 
low unless L is very small. The case of thin films (Fig. 
4b) is slightly more promising with larger efficiencies at 
longer lengths, but this requires small film thicknesses. 
Also note that the large AT over the small lengths L 
where reasonable efficiency is observed in Fig. 4a and 4b 
would lead to very large heat fluxes that further limit the 
usefulness of these geometries. Thus, for the bulk or thin 
film geometries to be useful, a much larger gT is needed 
so they can operate at longer lengths, and this would 
require other or new materials with better properties for 
SCL thermoelectrics. 

The situation is significantly more favorable if one con- 
siders nanowires, as shown in Fig. 4c. In that case, SCL 
currents are enhanced due to the scaling £ 2 ~ Rr 2 L~ l : 
in fact, gT becomes independent of length, and scales as 
R~ 2 , with R being naturally small for nanowires. This 
leads to a large efficiency which exceeds the 10% value for 
nanowires less than 20nm in radius. This efficiency could 
be even larger if other effects predicted and measured for 
nanowires were included; for example, from Eq. (1171) . gT 
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FIG. 4: Thermoelectric efficiency at T h = and AT = 

300A" for (a) bulk material as a function of length, (b) thin 
film material as a function of length for two thicknesses, and 
(c) nanowire as a function of radius. In (b), the curves are 
plotted for t < L/10 to satisfy the criterion t <C L (see Eq. 
[lil- 



dcpcnds sensitively on the bandgap, which increases with 
decreasing diameter. Similarly, the thermal conductiv- 



ity is reduced with decreasing nanowire diameter, which 
could also make gT larger. It should be noted that the 
electrical currents needed to achieve the efficiency plotted 
in Fig. 4c are less than a fiA per nanowire and require 
voltages on the order of a volt, which is easily achieved 
without damaging the nanowires or the contacts (l9j. 



IV. CONCLUSION 

In summary, we considered the efficiency of thermo- 
electric materials in the space-charge-limited regime and 
found that it depends on a single dimensionless param- 
eter gT, in analogy with conventional ohmic materials. 
When applied to bulk, thin film and nanowire geome- 
tries, we find that nanowires are the most promising to 
harness SCL transport. This work provides a new path 
for improving the performance of thermoelectric materi- 
als, and suggests the exploration of new thermoelectric 
materials with properties conducive to SCL transport. 
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